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Abstract 

The lattice model of principal chiral field interacting with the gauge fields, which 
have no kinetic term, is considered. This model can be regarded as a strong coupling 
limit of lattice gauge theory at finite temperature. The complete set of equations 
for collective field variables is derived in the large ./V limit and the phase structure 
of the model is studied. 
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1 Introduction 



The Kazakov-Migdal model Jl] was originally proposed as a model of induced QCD. 
Although this idea suffers from several problems @, there is another motivation for study 
of this model and it's generalizations - the usual for matrix models interpretation in 
terms of random surfaces, 2d gravity and noncritical string theory. The Kazakov-Migdal 
model significantly simplifies in the large N limit and can be treated by saddle point 
methods || £| , or by means of loop equations || . In this paper the modification of the 
Kazakov-Migdal model with the Hermitean matrix field replaced by the principal chiral 
one is considered: 



£ = Ntr 



±- (V^)t(V^) - V(g) - Vtf) 
2 1 



= d^g - [A^g\. 
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There g(x) is U(N) (or SU(N)) matrix, while A^(x) is u(N) (or su(N)) gauge field. 

This model, defined on a D-dimensional lattice, can be regarded as a strong coupling 
limit of the (D + l)-dimensional lattice gauge theory at finite temperature @, |7|, U- It 
can be shown in a most simple way within the Hamiltonian framework |J . In the strong 
coupling limit only the electric part of the Kogut-Susskind Hamiltonian M 



X,/l 



contributes to the partition function 



Z 



9 □ 



Tre T . 
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There ^ujx) is invariant Laplacian acting on the group variable U^x), attached to the 
link (x, x + fj) of a D-dimensional lattice. If it were not for the gauge invariance, the 
partition function ( |1.3|) would factorize in the strong coupling limit on the product over 
the links of the heat kernels on the group manifold, independently integrated over the 
gauge fields. However, one should first average over all gauges: 



UDU^x) UDg(x) UK[U ll (x),g\x)U fl (x)g(x + ^ 
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where K(g, h\r) is a solution to the heat equation 

^K(g,h\r) = ^A g K(g,h\r) 

with initial condition 

K(g,h\0) = S(g,h). 
An explicit expression for K(g,h\r) in terms of a character expansion is 

K(g,h\r) = ^ dim R e ~^Xr(9^), 

R 
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where the sum is taken over all irreducible representations of the gauge group, with dim/?, 
cr and xr being the dimension, quadratic Casimir and the character of the representation 
R. 

The partition function ( |1.4| ) is nothing that the lattice regularization of ( |1 . 1| ) with zero 
potential and 7 related to the gauge coupling g and the temperature T by 

9 2 N , \ 

7=V* (L8) 
This model, with exponential Boltzmann weight instead of the heat kernel one and a 
kinetic term for the gauge fields added, was studied both numerically and in a mean field 
approximation for D = 3,4 and N = 2,3 |TD[. The purpose of the present paper is to 



derive the large N equations of motion for collective variables in spirit of Refs. |Tl| , [12]1 and 
to study the phase structure of the model in the large N limit. 

It is well known that at some critical temperature lattice gauge theories undergo the 
deconfining phase transition || [13| associated with the center group symmetry breaking. 



In terms of the above model this is the symmetry with respect to the multiplication of 
all g(x) by an element of the center of the gauge group (this symmetry is spoiled by 
an addition of the potential term). At low temperatures (large 7) the model is in the 
symmetric (confining) phase and ^ tr g(x)"j=0, while at some critical value of 7 the 
symmetry is violated and the chiral field gains a nonzero expectation value. It will be 
shown below, that, at least for sufficiently large D, confining and deconfining phases are 
separated by an intermediate one. 

At high temperatures the fluctuations of the chiral field are suppressed, as 7 is small, 
so it can be expanded around unity 

g (x) = e^ {x) ~ l + 2 7 $(x) (1.9) 

and (|1.4j) reduces to the partition function for the Kazakov-Migdal model with critical 
potential M. 



2 Large N equations of motion 



Consider the partition function (|1.4j) with a potential term added: 

Z = [ Y[Dg(x)Y[DU,(x)e-^ Nt "l v ^ +v ^^ 

X X,fl 

l[K(g(x),U fl (x)g(x + fx)Ul(x) 
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7 



(2.1) 
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Integration over the gauge fields goes independently on each link, so, fixing the diagonal 
gauge: g(x) = diag( e mfc ( x )) and integrating out the link variables, one is left with an 
effective action for the eigenvalues of g(x), which contains only local interactions: 



eff 



E 



NY,U(a t (x))-J2\ogsm 

i i<j 



2 ai{x) — ctj [X ) 



2 



D 



log I (a(x),a(x + 



where 



U(9) = V(e ie ) + V(e 
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the second term in (2.2) comes from the decomposition of invariant measure: 



A? 



Dg = DnY[daiJ 2 (a), J(a) = ]Jsm 



Oti — OLj 
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(2.4) 
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for g = Qe ia Q\ and 



I(a,a'\r) = J DU K(e ia ,U e ia 'U ] \t). 



(2.5) 



It is convenient to introduce the density of the eigenvalues of the matrix g(x): 



N 



MM) = —^2S(9 - ai(x)) . 

iv i=l 



(2.6) 



In the large A" limit p(6, x) becomes smooth 27r-periodic function, normalized to unity on 
the interval (— n, n). For simplicity we restrict our consideration to the case of U (N) gauge 
group. For SU(N) the constraint J*^ d99p(9, x) — should be imposed on the eigenvalue 
density and the Lagrange multiplier ensuring this condition should be introduced. This 
is achieved by adding the term X(x)9 to the potential (|2.3j ). 

The classical equations of motion for (|2.2|) are exact in the large iV limit and can be 
written in terms of the eigenvalue density as follows: 



D 



U'(9) + (D-l)G(9,x)= Yl ^(M7), 



-D 



where 

and we have defined 
i^(M7) 



G(9,x) = p d9'p(9',x) cot 



1 d 



9-9' 



Fl^j l0g 1 {a{x) > a{x + A l7) + l G{e > X) - 
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The one-link integral (|2.5[) is the unitary analog of the Itzykson-Zuber integral and 



can be treated by the same method, proposed in Ref. [14] for finite iV and used to calculate 



the large iV limit of the Itzykson-Zuber integral in Ref. [|15|]. This method exploits the 
relation of the matrix quantum mechanics to that of the free nonrelativistic fermions |16| . 
Really, K(g,h\r) is an imaginary time transition function for the free Hamiltonian on 
the U(N) group manifold. After the calculations, similar to that of Ref. [II] , it can be 
shown, that J(a) J (a') I (a, a'|r) is the transition function for noninteracting fermions 



3 



on a circle. From which it immediately follows, that it can be written in a form of a Slater 
determinant: 



I(a, a'W) 



detfcj $ 



const 



2vr 



IT 

2ttN 



(2.10) 



J(a)J(a') 

where $(z, r) is the Riemann theta function. 

In the large N limit the fermions behave semiclassicaly, so the transition function for 
them is dominated by a single classical trajectory given, in terms of the density a(9, r), 
by Euler equations for an ideal fluid [n, I2|: 



do- d . . 

TTr + Te [as) = °' 



<9s 9s 



da 



tt" + s— - ix o— = 



dr 
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(2.11) 



(2.12) 



The minus sign before the last term in ( |2.12|) is because the time is imaginary. The 
function s(9,t) is a velocity of the Fermi fluid. Eqs. Q2.11|) , (|2.12|) , of course, can be 



derived directly by a reduction of the heat equation for (|2.5| ) to the Hamilton-Jacoby 
equation for the collective field Hamiltonian by the calculations, analogous to that of Ref. 

To find ( p.9| ) one should solve ( |2.11| ), ( |2.12| ) with boundary conditions 



Then for (|2.9|) one has: 



(2.13) 
(2.14) 

(2.15) 

This formula follows from the definition ( |2.9|) , because the variation of the classical action 
with respect to the initial condition for the density Q2.13j ) is equal to it's conjugate variable, 
which is a potential of the velocity fT7|| . Reversing time in ( 2.11| ), ( 2.12 ) one also finds: 



a(9, 0; x, li) = p(9, x), 
a(9, 7; x, fi) = p(9, x + fi). 

F fl {9,x\j) = s{9,0;x, f x). 



F_ M (0,x|7) = -s(9,r,x - ix, n). 



(2.16) 



Thus in the large N limit the model reduces to the set of differential equations ( ^.11| ), 
( ^.12| ) to be solved on each link with boundary conditions for the density given by ( ^.13| ), 
fl2.14j ) and the sum of the jumps of the velocities at each lattice cite given by the l.h.s. 
of eq. ( |2.7|) . Eqs. ( [2.11D , Q2.12 ) are integrable, although in an implicit form, and the 
problem is, in principle, reducible to the purely algebraic one jl5| , differing from that for 
Hermitean matrix model in the periodic boundary conditions for 9 and the kernel in the 
integral (|2.8|) (the Cauchy one versus the Hilbert kernel in the Hermitean case). 



In one-dimensional case the gauge field in ( |1.1| ) can be absorbed by a gauge transfor- 
mation, so D = 1 model is equivalent to the unitary matrix quantum mechanics, solved 
in the large iV limit in Refs. |12|, 113]. It is instructive to reproduce these results from the 
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above lattice equations. One can suppose that r in eqs. Q2.11|) , (|2.12|) changes from — oo 
to +00, while at the points t = nj the velocity has a jump according to ( |2.7|) , ( |2 ■ 1 5|) , 
( p,16| ). However, in the continuum limit 7 and U\9) scales as 07 and all' (9), respectively 
(a is the lattice spacing), so the discontinuities in s(9,t) can be absorbed by adding the 
term W{9) to the r.h.s. of eq. ( gl2|) . The magnitude of a deviation of the solution 



of resulting equations from the exact one on the time interval (ncry, (n + 1)07) is of the 
order (07) 2 and can be neglected as a — > 0. Thus one obtains in the continuum limit: 

| + 7|(H = 0, (2.17) 
dv dv 9 dp < 

where t is the real time, measured in an ordinary units (while r is measured in the units 
of 7). These are nothing that Jevicki-Sakita collective field theory equations ||12|| . 

In the multidimensional case the second term in the l.h.s. of eq. (|2.7|) do not vanish 
and is of order unity as a — > 0, so to obtain the continuum limit in the same way as in one 
dimension, the solutions should be considered with the velocity s(9, r; x, p) being by itself 
of order a -1 and 7 scaling as a 2r y. Then the arguments, similar to that used to obtain 
( p.lTp . ( p. 18 ), give the following equations in the continuum limit (in the Minkowski 
space): 

3i}^ 3i]^ 

^ + ^ + (D-l)G + U> = . (2.20) 

As in the case of Hermitean matrix model, it can be shown, that translationally invariant 
solution to these equations is unstable the spectrum containing an infinite number 
of tachyons. Thus one has to look for other possibilities to take the continuum limit. In 
this respect the study of the phase structure of the model is of fundamental importance. 



3 Strong coupling phase 



Consider the model ( |2.1| ) without a potential term. Equations of the previous section 
always have a trivial solution p(9, x) = ^- in this case (with s = in ( [2. lip , ( |2.12| )). This 



solution is valid at large 7. However, the instability in the continuum limit at 7 = 
shows, that at some 7 C the strong coupling solution becomes a maximum, rather than a 
minimum, of the effective action. 

To find the spectrum of excitations about the strong coupling solution we write 

p{e,x) = j- + j-^a n (x)e infi , a* n = a. n (3.1) 

-/ n^O 

cr(9,T;x,n) = — + — a n (r;x,p)e me , a* n = a_ n (3.2) 
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(3.3) 



s(9,r;x,p) = £ (3 n (r; x, p) e in0 , j3* n = (3_ n 

and linearize the equations of motion in a n , a n and j3 n . The solution to eqs. Q2.11| ), (|2.12 
then reads: 

a n (r;x,p) = a^(x,p)e^ + a~(x, p) e (3.4) 
% 
2 



f3 n (r;x,p) = - a^(x,fj,)e2 — a n (x,p)e 2 



(3.5) 



Substituting these expressions in the boundary conditions ( [2.13| ), ( |2.14] ), ( |2.15| ), ( |2.16| ) 
and ( |2.7p one obtains, after a simple algebra, the equation for the Fourier coefficients of 
the eigenvalue density: 



D 

E 



a n (x + /i) — 2 cosh 



Try D — 1 n.7 

sinh — 

£> 2 



a n (x) + a n (x - /i) 



0. 



which leads to the following mass spectrum: 



Ml = D cosh ^ - (D - 1) sinh 

where = X^ =1 cosh for the solution of ( PTB{ ) of the form exp (Z)^=i A^x M 
At the critical point 

7c = 21og(2D-l) 



(3.6) 



(3.7) 



(3i 



the lowest excitation becomes massless, and the strong coupling solution becomes unsta- 
ble. In principle this critical point can be used to construct the continuum limit. However, 
it describes only one complex scalar particle with the mass 



m 



[D-l 



7 ~7c 

n2 ' 



(3.9) 



while the masses of other excitations remain of the cut-off order. It is worth mentioning, 
that an obstacle to take the continuum limit at 7 C may occur due to the possibility for 
the model to undergo a first order phase transition into the weak coupling phase before 
the point (|3.8|) is reached (see the next section). 



4 Weak coupling phase 



As 7 — > the expansion Ql.9|) becomes legitimate, so one expects that the eigenvalue 
density p(6, x) is peaked about zero. Thus, in the zero approximation, the cotangent 
in ( |2.8| ) can be replaced by and the equations of motion become the same as 

for the Kazakov-Migdal model with the action containing only the derivative squared 
term. It is known, that such model has a stable translationally invariant solution [J3j]. The 
fluctuations about it contains one massless mode |2(J related to the shifts of the center of 



the eigenvalue density. It is straightforward to show (in the general case) that, if p(9,x) 
solves the equations of Sec. 2, then p[Q + #0(3;), x) also does, provided that 9q(x) satisfies 



6 



lattice Laplace equation. So the Goldstone mode simply decouples and can be disregarded 
(one may exclude it considering SU(N) model). 



Now let us take into account the next term in the expansion of the cotangent in ( |2.8| ) 
(the eigenvalue density is assumed to be an even function of 9): 

p J d9'p(9', x) cot 6 -^L = 2 p J d9'f^- - l -9 - . . . (4.1) 



After the substitution of fl4.1|) in the l.h.s. of eq. (|2.7|) the equations of motion become 



equivalent to that for the Kazakov-Migdal model with a mass term added. Note, that 



this mass is negative. Moreover, the residual term in ([4.1|) is given by the convergent, 
as \9 — 9'\ < 2tt, series in the odd powers of 9 — 9' with negative coefficients. Thus an 
account of the corrections coming from the nonlinearity of the chiral field leads to the 
p-dependent upside-down potential and ( |4.1|) gives an upper bound for it. For sufficiently 
weak potential an interaction with the gauge fields stabilizes the fluctuations and provides 
a minimum for an effective action, but at some critical 7* the effective potential becomes 
strong enough to make a weak coupling solution unstable. An upper estimate on 7* can be 
obtained considering the Kazakov-Migdal model with the quadratic potential, an effective 
mass being given by (in the notations of Ref. @]): 

m 2 eff = 2D-±(D-l) % (4.2) 
the solution to which is known M: 



p{0)^\^-\^e\ (4.3) 



m 2 eff (D-l) + DJmt ff -4(2D 



This solution is stable until roL = 2\/2D — 1, when (|4.4| ) becomes complex and the 



massless excitation appears in the spectrum pOfl . This gives the following bound on 7*: 



12(D - \J2D - 1) 
7 * < D _i • ( 4 - 5 ) 

It is worth mentioning that for D = 2, 3 the width of the support of eigenvalue density 
( |4.3| ) becomes larger than 2tt for smaller value of the coupling. 

There are, in principle, three possibilities for the critical behavior: 



i) 7* > 7 C . In this case strong and weak coupling solutions can coexist and the first order 

phase transition takes place at some 7 between 7 C and 7*. 

ii) 7* = 7 C . The phase transition is of the second order. 

iii) 7* < 7 C . This means that an intermediate phase exists, separated from strong and 

weak coupling ones by the second order phase transitions. 
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The third possibility is realized in the large D limit. Really, 7 C grows as 21ogD, while 
7* tends to 12 as D — > oo (the estimate (|4.5| ) gives an exact large D asymptotics of 
7*, because the width of the eigenvalue distribution ( }4.3| ) behaves at the critical point 

as 4 and corrections to ( |4.1| ) vanish). However, using an upper bound ( |4.5| ) it 

is possible to advocate the existence of an intermediate phase only for D > 86. It is 
worth mentioning that the value of 8, at which the eigenvalue density (|4.3|) with critical 
miff and 7, given by an upper bound ([15]), turns to zero, is, nevertheless, not small - 



6max — 0.66-7T for D = 87, so the approximation ( |4.1| ) remains rough at the point of the 
phase transition even for such large D. Thus more accurate treatment is necessary to 
distinguish between the possibilities i)-iii) for low dimensions. 

For D = 3, 4 the first possibility is favored by the Monte Carlo simulations of Ref. [|K]] 
and the large amount of numerical data on the finite temperature gauge theories. However, 
the second possibility do not contradict the numerical results, as we are dealing with the 
boundary of the phase diagram in g 2 , T, and the line of the first order transitions may 
terminate at the point of the second order one. The third possibility is also not excluded, 
because an intermediate phase can disappear at some finite, but large N, and thus be 
invisible in the Monte Carlo simulations for N = 2, 3. 



5 Conclusions 

If we discard the possibility of the first order phase transition between the strong and 
the weak coupling phases, which is legitimate, at least for large D, as discussed in Sec. 
4, the continuum limit of the model can be constructed. Unfortunately, corresponding 
continuum theory contains only finite number of degrees of freedom, if we approach the 
critical point from strong or from weak coupling phase, and thus cannot describe extended 
objects, which was one of the motivations for study of the model. However, this may not 
be the case for the model with a potential term added, or if one takes the continuum limit 
in the intermediate phase. It is also interesting to understand the origin of this phase from 
the point of view of thermal lattice gauge theory. To clarify these questions it is desirable 
to know an exact solutions to the model in the weak and the intermediate phases. To 
this end, it would be interesting to develop an alternative approach, based on Schwinger- 
Dyson equations, which proved to be useful in the context of the Kazakov-Migdal model 
§• 
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Note added 

When this work was being prepared for publication, there appeared a paper by D.V. Boula- 



tov 21 , which also includes the treatment of the gauged chiral field model at large N. 
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